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4 Microgravity propellant dynamics continue to offer a formidable 
modeling challenge for the computational fluid mechanics community 

•Analytical approaches to prediction of bulk behavior, e.g. tank forces 
and moments, degrades in accuracy below Bo-IO; small perturbations 
only 

• Flows dominated by surface tension; curved interfaces; fluid-wall 
contact angles approximately zero 

• Many semi-analytical or empirical methods that correlate well with 
theory rely on quasi-steady-state parameters and cannot accurately 
predict effects of transient flows, e.g. throttling, thruster pulses 

• Momentum transfer due to fluid must be computed accurately for 
simulation-based verification 


4 Our research explores the applicability of the iattice Boitzmann 
method (LBM) to modeiing of cryogenic propeliant dynamics in 
microgravity 
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The Lattice Boltzmann Method 
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4 The lattice Boltzmann method (LBM) is an emerging approach to CFM using an 
explicit temporal and spatial discretization of the continuous Boltzmann equation: 

^ + a^V^/ = f2(/) 

• Describes evolution of particle distribution functions, e.g. density distribution 

• Regular Cartesian discretization of a 2XD position-velocity phase space on a lattice 

• In discrete form, the lattice Boltzmann equation is given by 


f(x/, + e5t,to + 5t) = f(x/,,to) - A(f(xjfc,to) - r"^(xfc, u^, p(xfc))) 


Adjacent Next time Current Relaxation 
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p(xfc)ufc = y^e^/^(xfc), 
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Numerical Implementation and 

Limitations 
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♦ The fluid dynamics are propagated in two steps 
over multiple layers 

1 . Advection (copy fluid density to adjacent cell: 
“streaming”) 

2. Collision (simulate collisions by relaxing 
toward Maxwell equilibrium: “relaxation”) 

♦ Proper choice of units (time=space=1) and 
periodic lattice: no actual data copy 

• The streaming step is done using pointer 
arithmetic: fast 

• Data locality (only need knowledge of 
adjacent cells): parallelization 

4 Limitations: 

• Constant (wall) temperature: isothermal flow 

• Effective speed of sound is related to the 
lattice physical parameters 

- Compressibility error increases as M>0A 

- Incompressible flow approximation is valid 
only for steady flow at low velocities 

• Stability decreases as timestep decreases at 
a fixed kinematic viscosity 
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Boundary Conditions & Body Forces £ 





4 Wall boundary conditions are straightforward to implement 

• So-called “bounce back rule”: invert velocity distributions at boundary 

• Reconstruct unknown distributions by storing in boundary for one timestep 

• Not a hydrodynamically accurate BC, but simulates no-slip wall 

• More accurate BCs are required for free slip, walls with high curvature, inlet flow, etc. 


♦ Relaxation operation and body forces 

• Body forces implemented using Kupershtokh exact difference method (EDM) 

• Shift equilibrium distribution under action of force such that lattice remains in equilibrium 

• Relaxation operation uses multiple relaxation time (MRT) scheme accounting for 
variation in kinematic viscosity with optional subgrid turbulence model 


t — St t 



Half-Way Bounce 
Back BC 
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Instrumented 6 P at depth = 0.494 m 
1-cos() acceleration profile to 0.100 g, max 6 1 = 0.207 ms 



simulation time (sec) 
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Flow & Wall Wetting 
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4 LBE models allow incorporation of muitipie phases uniformly in the iattice 

• Phase separation explicitly depends on temperature and a real gas EoS 

• We use the Carnahan-Starling EoS corrected for the target conditions (LO 2 @ 94K) in a 
Shan-Chen like pseudopotential model 


p = pRT 


1 + 6p/4 + ihp/A)^ — {hp/A)' 

(1 - hp/Af 



• Phase segregation approximately obeys Maxwell construction (“mechanical stability”) 


4 Parametric wail wetting model allows tuning of free surface contact angie 
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Simulating Droplet Dynamics 
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4 The behavior of droplets in microgravity is dominated by surface tension 

• Droplet dynamics provide a useful verification case due to analytical solutions 

• Oscillation of the free surface can be predicted by Lamb’s equation 

• Effective surface tension can be determined from frequency 
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*Frohn, A., and Roth, N. Dynamics of Droplets. Springer, 2000 
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Propellant Slosh (Low Bo) 


4 Propellant sloshing dynamics are of fundamental concern in spacecraft dynamics 

•Analytical solutions exist for axisymmetric containers at high Bo 

• Below So=1000 modified analytical models must be used 

- Stable flow only (continuous free surface) and small perturbations 

• CFD solutions required for turbulent flow, transient phenomena, PMD simulation, etc. 



4 LBE method verified using a fiow regime near anaiytical limit 

• 0.15 m cylindrical LO 2 tank with hemispherical ends @ 0.001 g (Bo=20) 

• Lattice size = 346^ (4.1 MB) 

• Free decay from initial condition with acceleration 15 degrees from symmetry axis 

• First mode frequency matches analytical predictions very well (/=0.055 Hz) 
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Fluid Slosh (Microgravity) 
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4 LBM approach used for small domain simulation in microgravity 

• 0.1m spherical LO 2 container in microgravity 

• High-g initial condition (settled fluid) and Og at t=0 (20% fill fraction) 

• Random perturbing acceleration (-3 dB @ 0.2 Hz, gp^g^O. 00003) 

• Lattice size = 266^ (2.4 MB) 


density <kgAn3) (mean (-178.1. 20 6-1096 8) 
ln>e ■ 8 357 sec. dt-16 614 ms 



dens4y (kgAn3) (mean (-178.1. 14 5-1 110 4) 
time - 20 003 sec. dt-16 614 ms 



density tkgAnS) (mean (-178.1. 13 2-1 1 18 1) 
t»ne ■ 60 010 sec. dt-16 61 4 ms 
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density (kgAn3) (mean (-141.3, 10.3-1093.0) 
bme - 90 010 sec. dt-15 109 ms 


density (kgAnS) (mean (-141 3. 11.9-1103 4) 
tune - 130 004 sec dt-15 965 ms 


density (kgAn3) (mean (-141 3. 14.0-1097.0) 
txne- 150 008 sec. dt-16 614 ms 


denuy (kg/lT)3Mmean (-141.3.112-1117 8) 
time- 180 013sec. dt-l66i4ms 
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4 Some promising resuits have been obtained in the application of the LBE to 
propeliant sloshing dynamics and reiated microgravity fiuid phenomena 

•The LBE method has promise for simulation of thermodynamically consistent multiphase 
flows in cryogenic propellants 

• Fundamental stability limitations remain for very low kinematic viscosities (real fluids!) 

• Recent progress includes stability enhancement via adaptive time-stepping and 
improvement in surface tension model using multirange pseudopotential 

- Overcomes some stability and surface tension limitations by improving isotropy 

4 Production simulation of cryogenic flows will require incorporation of thermal effects 

• Thermal LBE codes are emerging and show some promise 

• Important to capture convective phenomena, for example, for modeling of long-term fluid 
circulation in propellant storage systems 

4 Ongoing work extends the present proof-of-concept model to a 3D code 

• Parallelization opportunities may allow simulation of CFD-in-the-loop with spacecraft 
GN&C 6-DoF simulations, for example, using GPU computing 

• Possible validation opportunities using existing microgravity fluid experiment data 
collected aboard International Space Station (ISS) 

• 3D code targeted for NASA Exploration Upper Stage (EUS) ullage settling and propellant 
management studies 
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